Electron transport in a Pt-CO-Pt nanocontact: First-principles calculations 
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We have performed first-principles calculations for the mechanic and electric properties of pure Pt 
nano-contacts and a Pt contact with a single CO molecule adsorbed. For the pure Pt contacts we see 
a clear difference between point contacts and short chains in good agreement with experiments. We 
identify a tilted bridge configuration for the Pt-CO-Pt contact, which is stable and has a conductance 
close to 0.5Go (Go = 2e 2 /h), and we propose that this structure is responsible for an observed peak 
at 0.5Go in the conductance histogram for Pt exposed to a CO gas. We explain the main features 
of the transmission function for the Pt-CO-Pt contact, and show that the conductance is largely 
determined by the local d-band at the Pt apex atoms. 

PACS numbers: 73.63.Rt, 73.20.Hb, 73.40.Gk 



I. INTRODUCTION 

Advances in the experimental techniques for manip- 
ulating and contacting individual atoms and molecules, 
and the vision of using simple organic molecules as the 
basic building blocks in electronic devices, have recently 
intensified the interest for electron transport in nano- 
scale contacts i2i2i4 On the theoretical side a number of 
first principles methods to describe the electrical prop- 
erties of realistic atomic-sized junctions have been devel- 
oped .SjejTjSjQ Marry of these methods combine a single- 
particle description of the electronic structure extracted 
from density functional theory (DFT)^Si with a non- 
equilibrium Green's function formalismAl to calculate the 
current. The various schemes mainly differ in the choices 
of basis sets and the way in which the coupling between 
the molecule and macroscopic leads is taken into account. 

Regardless of the details of the implementation, the 
LDA/GGA DFT based transport schemes show the same 
general trend: the calculated conductance for molecu- 
lar contacts in the low-conducting regime, i.e. contacts 
with a conductance much lower than the conductance 
unit G = 2e 2 /h, is 2-3 orders of magnitude higher than 
corresponding experimental values A number 
of reasons for this discrepancy have been proposed, in- 
cluding correlation effects not captured by the mean-field 
approachiii& and differences between the atomic struc- 
tures used in the calculations and those realized under 
experimental conditions piiiSi On the other hand, a most 
satisfactory agreement between DFT based calculations 
and experiments is found for the conductance of homo- 
geneous metallic point contacts and monatomic wires. 
For these systems, both the size of the conductance, 
the number of conductance channels, and the so-called 
conductance oscillations are well reproduced by calcu- 
lations jSiiSLSLSiiSSiSiLSiiSiSSi The same good agreement 
between theory and experiment has been found for a het- 
erogeneous molecular junction consisting of a single hy- 
drogen molecule captured between Pt electrodes iSLSSiSS 
Common to all these systems, for which DFT provides 
a description of the transport in agreement with experi- 



ment, is that the coupling to the metallic leads is strong 
and the conductance is high, i.e. on the order of IGq. 

There remains much to be learned about the reasons 
for these apparent trends in DFT based transport calcu- 
lations and it is therefore desirable to extent the knowl- 
edge of systems, for which transport properties can be 
treated within the framework of DFT. 

In this paper we present DFT calculations for pure 
Pt contacts and a Pt-CO-Pt system, which according to 
experiments have a conductance on the order of IGq and 
thus belongs to the group of systems for which previous 
DFT transport calculations have been successful. We 
have calculated the total energy and conductances of pure 
Pt contacts and chains as well as for Pt contacts with a 
single CO molecule adsorbed. Our main result is the 
identification of a certain "tilted bridge" configuration 
for the Pt-CO-Pt contact which is energetically stable 
and has a conductance close to 0.5Go in agreement with 
recent experimental results^ We find that the transport 
properties of the Pt-CO-Pt junction to a large extend are 
determined by the properties of the bare Pt electrodes. 
For this reason we put some emphasis on verifying the 
ability of our method to reproduce key characteristics of 
the transport properties of pure Pt contacts and chains. 
Also here we find good agreement with experiments. 

Mechanically controlled break junction experiments 
performed at cryogenic temperature on pure Pt samples 
show, that as a Pt contact is pulled apart a structure with 
a characteristic conductance of around 1.5Go is formed 
in the last stages before the contact breaks. This is in- 
ferred from conductance histograms^2*2iiS which show 
a pronounced peak at this value. In addition to the 
peak at 1.5Go, many histograms on Pt contain a smaller 
and broader peak at around 2.1 Go. The two peaks are 
believed to correspond to chains and atomic point con- 
tacts, respectively. The fact that the peak at 1.5Go is 
higher than the peak at 2.1Go is explained by the sup- 
pression of point contacts by the formation of chains. 
Experimental evidence for this hypothesis comes from 
conductance histograms recorded as the broken contacts 
are brought back into contact - the so-called return his- 
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tograms. Such histograms contain no contributions from 
chains, and show a single peak at 2.1G . 31 Our calcula- 
tions predict that Pt point contacts have a conductance 
of (2.0 — 2.3)Go, whereas short Pt chains have a con- 
ductance of (1.3 — 2.0)Go in good agreement with the 
experimental findings. 

When the Pt contact is exposed to a CO gas, the peaks 
characteristic for pure Pt disappears from the conduc- 
tance histogram and are instead replaced by two peaks at 
~ 0.5Go and ~ I.IGq. The physical origins of these peaks 
have not been identified experimentally Recent measure- 
ments on magnetic as well as non-magnetic metal point 
contacts^., showed a fractional conductance of 0.5Go, 
and this was interpreted as the lifting of a spin-degenerate 
conductance channel. However, none of these results 
could be reproduced in experiments by C. Untiedt and 
co-workers. 30 Instead, they suggested that the reported 
fractional conductances could result from CO contami- 
nation. Regardless of whether CO is the source of the re- 
ported fractional conductances, the question of the phys- 
ical mechanism of the peak at 0.5Go in the Pt-CO-Pt 
histogram still remains: is it a spin effect or does it have 
some other origin ? 

On the basis of our calculations we propose that the 
0.5Go peak in the Pt-CO-Pt histogram is due to a con- 
figuration where a single CO molecule provides a "tilted 
bridge" between two Pt apex atoms. We find that this 
structure has a conductance just below 0.5Go for a wide 
range of electrode displacements leading to a plateau in 
the conductance trace. By carrying out a Wannier func- 
tion analysis we identify the current carrying states of the 
molecule as the 2ir* CO orbitals. However, the main fea- 
tures of the transmission function are not determined by 
these orbitals but rather by the Pt d-band at the apex 
atoms. The fact that the Pt-C bond is much stronger 
than the Pt-0 bond allows us to simplify a resonant level 
model for asymmetric coupling, and obtain a simple de- 
scription of the transmission function in terms of a single 
d-like Pt orbital, the energy level of the 2ir* CO orbital, 
and a coupling strength. All the parameters of the model 
are extracted from the first principles calculations. 

The paper is organized as follows. In Sec. [H] we out- 
line our Wannier function based transport scheme. In 
Sec. IIIII we present the results for simulated conductance 
traces of pure Pt point contacts and short Pt chains. In 
Sec. IIVI the Pt-CO-Pt system is investigated by simu- 
lating a conductance trace and identifying the origin of 
the main features of the transmission function. Finally, 
Sec. contains a summary. 



II. METHOD 

In this section we briefly review the computational 
methods used in the present study. All total energy cal- 
culations are performed using a plane wave implementa- 
tion of density functional theory^* The nuclei and core 
electrons are described by ultra-soft pseudopotentials 35 , 



and exchange and correlation is treated at the GGA level 
using the PW91 energy-functional^ The Kohn-Sham 
(KS) eigenstates are expanded in plane waves with a ki- 
netic energy less than 25Ry. Optimizations of all the 
considered structures have been performed, until the to- 
tal residual force is below 0.05eV/A. 

The electrical conductance is evaluated within the 
Landauer-Butikker formalism. Thus the system is di- 
vided into three regions: a left lead (L), a right lead (R), 
and a central region (G). The leads are assumed to be 
perfect conductors such that all scattering takes place 
in G. The linear response conductance due to scatter- 
ing upon G is given by G = GqT(sf), where T{ep) is 
the elastic transmission function evaluated at the Fermi 
energy. Here we evaluate T{e) from the self-consistent 
KS Hamiltonian, which we represent in terms of a basis 
consisting of maximally localized, partly occupied Wan- 
nier functions (WFs). The WFs are constructed accord- 
ing to a recently developed method^LiSiSi in such a way, 
that the KS eigenstates are exactly reproducible in terms 
of the WFs within some specified energy window. The 
energy window is typically selected to include all eigen- 
states up to around 4eV above the Fermi level, and thus 
the plane wave accuracy of the original DFT calculation 
is retained in the subsequent transport calculation. The 
use of a localized basis set allows us to calculate T(e) 
using the general non-equilibrium Green's function for- 
malism. 40 This formalism is formulated in terms of or- 
thogonal basis functions, but was recently generalized to 
the case of non-orthogonal basis functions^i 

Furthermore, the formalism is generally valid for a fi- 
nite bias voltage, however, in this study we shall focus 
on the linear response conductance. 



III. PT POINT CONTACTS AND CHAINS 

In this section we investigate the electrical properties 
of pure Pt contacts and short Pt chains between bulk 
Pt electrodes. This study serves a dual purpose. First, 
it provides a theoretical justification for assigning the 
two peaks at 1.5Go and 2.1Go in the conductance his- 
tograms for Pt to chains and point contacts, respectively. 
Secondly, it allows us to test the ability of the calcula- 
tional scheme against well established experimental re- 
sults and other computer simulations before applying it 
in the study of the Pt-CO-Pt contact. 

The Pt contacts are modeled using supercells by two 4- 
atom pyramids oriented towards each other and attached 
to Pt(lll) surfaces containing 3x3 atoms in the surface 
plane. Two different structures are considered: (i) A 
point contact, where the apex atoms of the pyramids are 
in direct contact, (ii) A 1-atom chain, where a single 
Pt atom is inserted between the apex atoms of the pyra- 
mids, see Fig. ^ hi order to ensure that the effective 
KS potential has converged to its bulk values at the end 
planes of the supercell we include 3-4 atomic layers on 
either side of the pyramids. 
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FIG. 1: Supercells used to model the two considered struc- 
tures: (i) A point-contact, and (ii) a 1-atom chain. The elec- 
trode displacement d z is defined as the distance between the 
(111) surfaces. 



bond between the two apex atoms is broken, which can be 
seen directly from the insets. The breaking of the central 
bond marks the onset of a structural relaxation, which 
affects both the conductance and force significantly. Be- 
yond this point the contact enters the tunneling regime 
and the conductance decreases exponentially. 

The simulated conductance trace is in good agreement 
with the experimental return histograms for Pt^ 1 which 
show a peak around 2.1Gn. Moreover, both the plateau 
around 2Gn as well as the rate of the exponential decay in 
the tunneling re gim e compare well with the calculations 
reported in Ref. |2(| 



For the total energy calculations, the Brillouin zone 
(BZ) is sampled by a 4 x 4 x 1 Monkhort-Pack grid4£ 
The transmission function is sampled over a 4 x 4 k-point 
grid in the two dimensional BZ of the surface plane. This 
sampling is crucial in order to avoid unphysical features 
in the transmission functions due to Van Hove singulari- 
ties associated with the quasi-one dimensional leads. 43 



A. Pt point contact 

By increasing the electrode displacement d z , defined as 
the distance between the fixed (111) surfaces (see Fig.0, 
and relaxing the pyramids at each step before calculating 
the conductance, we simulate the process of creating a 
conductance trace. The result for the point contact is 
shown in Fig.0 where the triangles denote conductances 
and the circles denote the total energies measured relative 
to the first configuration (d z — 10. 9A). 
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FIG. 2: Conductance (triangles) and total energy (circles) for 
the Pt point contact as a function of electrode displacement 
d z . 

For the first three configurations, the conductance 
stays just above 2Gn. The reason for the relatively 
weak changes in the conductance in this region, is that 
the elongation occurs quite uniformly over the pyramids, 
such that the distance between the apex atoms changes 
only slightly. At configuration four (d z — 12.07A) the 



B. Short Pt chain 

We have simulated the breaking of a one-atom Pt chain 
following the same procedure as for the point contacts 
described in the preceding subsection. The results are 
shown in Fig. [3] The calculated conductance trace has 
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FIG. 3: Conductance (triangles) and total energy (circles) as 
a function of electrode displacement, d z . 

a plateau at around 1.4Gn for small d z where the chain 
zigzags. As the contact is stretched further the conduc- 
tance rises to a conductance just below 2Gn, before the 
contact breaks at d z = 14.5A. At this point the struc- 
ture relaxes towards the surfaces, and the conductance 
starts to decay exponentially as the tunneling regime is 
entered. The correlation between structural relaxations 
and sharp changes in the conductance is a characteristic 
feature of the contact formation process. The effect has 
been observed experimentally for gold chains by mea- 
suring the conductances and forces simultaneously^ 4 ^ 
Furthermore, we notice that the increase in the conduc- 
tance just before the contact breaks, seen in both Fig. |3 
and Fig. is also observed experimentally, and is in fact 
characteristic of conductance traces for Pt and Al con- 
tacts. 32 In the case of Pt this behavior could be related to 
the linearizion of the chains which activates more conduc- 
tance channels^ A similar explanation has been given 
in the case of Al contactsiSiSi 

Finally, we stress that in order to obtain plateaus in a 
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simulated conductance trace and thus predict the occur- 
rence of peaks in a conductance histogram, it is necessary 
to allow the central atoms of the contact to relax in the 
elongation process. 



IV. PT-CO-PT CONTACT 

As discussed in Sec.[|] the controlled exposure of the Pt 
contact to a CO gas changes the conductance histogram 
completely: the peaks at 1.5Go and 2.1Go characteris- 
tic of pure Pt are replaced by peaks at 0.5Go and at 
1.1 Go- In order to understand the physical origin of the 
new peaks, we have carried out total energy and conduc- 
tance calculations using the same setup as for the pure 
Pt contacts discussed in the preceding section. 

The results are summarized in Fig. where we show 
the conductance (triangles) and total energy (circles) as 
a function of the electrode displacement d z . For the ini- 
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FIG. 4: Conductance (triangles) and total energy (circles) 
as a function of electrode displacement d z . The conductance 
trace is divided into an upper and a lower part, corresponding 
to the CO in the upright bridge- and tilted bridge configu- 
ration, respectively. A and B labels the configurations just 
before and just after CO has tilted, respectively. 

tial contact geometry and small d z , CO is most stable 
at the central Pt-Pt bridge of the contact, bonding sym- 
metrically with C to the two Pt atoms in an upright 
bridge configuration. It is seen that the conductance 
decreases linearly from 2.1Go to 1.5Grj as the upright 
bridge configuration is being stretched. At d z = 12. 8 A 
one of the C-Pt bonds breaks, and the contact relaxes 
to a tilted bridge configuration. This qualitative change 
in the atomic structure is clearly seen in the conduc- 
tance which jumps abruptly from ~ 1.5Grj to a plateau 
at ~ 0.5Go which then extends over almost lA. 

In order for the upright and tilted bridge configura- 
tions to be observed with a reasonable probability in the 
experiment, the binding energy of CO in these positions 
should be larger than that of CO adsorbed on the nearby 
sites, i.e. on the (lll)-surface sites. This is indeed the 
case, as the binding energy for the upright and tilted 



bridge is around -2.6eV, while that of CO on Pt (111) 
varies between -0.3eV and — 1.8eV, depending on the cov- 

4fi 47 

erageifiaai 

Whether the conductance plateau at 0.5Go in Fig. 0] 
will contribute to a peak in a conductance histogram, de- 
pends crucially on the stiffness k s of the Pt electrodes^ 
Each Pt electrode can be viewed as a spring connected 
to the CO contact, and in our setup only the relaxation 
of the atoms near to the CO molecule have been taken 
into account. These local relaxations are very important 
and responsible for the abrupt drop in the conductance 
curve. However, to obtain a more realistic conductance 
trace, the finite spring constant of the remaining elec- 
trodes should be taken into account, and this will de- 
form the cc-axis in Fig. 0] We make the assumption that 
at some points far away from the contact, the position of 
the electrodes are controlled and denote the separation 
between these points by L. The relation between L and 
the distance between the Pt surfaces, d z , can be obtained 
by solving the equation 



d z ) = 



dE{d z 
dd z 



(1) 



which expresses the force balance between the springs 
(left hand side) and the contact region (right hand side). 
The factor 1/2 is due to the fact that each springs are 
acting on both sides of the contact, with a stiffness of k s . 
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FIG. 5: Conductance traces for different values of the elec- 
trode stiffness, k s . The size of k 3 determines the length of the 
conductance plateau at 0.5Gn: a stiff electrode (large k 3 ) re- 
sults in a long plateau, while a soft electrode (small k s ) results 
in a short or no plateau. 



In Fig. El w e show the corrected conductance traces 
for k s between 1.8 eV/A 2 and 9.0 eV/A 2 . For k s be- 
low ~ 2 eV/A , the 0.5Go conductance plateau is com- 
pletely absent and the tilted bridge configuration would 
not contribute to the conductance histogram. However, 
when k s > 3 eV/A the 0.5Go plateau becomes clearly 
visible, and would lead to a peak at 0.5Go in the con- 
ductance histogram. Experimental estimates for k s do 
not exist for Pt electrodes, however, in the case of gold 
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values in the range 0.2eV/A to 4eV/A have recently 
been reported, on the basis of non-exponential distance 
dependence of the tunneling current. 49 

We therefore propose, that the observed peak in the 
conductance histogram for CO in Pt nanocontacts is due 
to the tilted bridge configuration. All our calculations are 
averaged over spin, and the special conductance of 0.5Go 
is therefore not related to spin. This is further illustrated 
below, where we study the transport mechanism in more 
detail. 

Finally, we mention that the peak around l.lGo in the 
experimental histograms, cannot be explained by our cal- 
culations for a single CO molecule in the Pt contact. We 
have also calculated the conductance of the CO bridge 
in the one atom chain in search of a structure that can 
explain the peak at l.lGo, however, we found conduc- 
tances quite similar to those obtained for CO in the Pt 
point contact. 



A. Conduction mechanism 

In this section we address the question of the physical 
origin of the fractional conductance of 0.5Go found for 
the tilted bridge configuration. Fig. [fj] shows a typical 
transmission function for the upright bridge configura- 
tion (dashed line) and tilted bridge configuration (full 
line) immediately before and after the bridge tilts - con- 
figurations A and B in Fig.Q] The dotted line denotes the 
transmission function obtained for configuration B with a 
localized atomic-like basis set taken from the DFT code 
Siesta^ The transmission function has been averaged 
over the same k-points and calculated using the same 
electron transport code as the WFs transmission func- 
tions. The basis set used is the Siesta "default" double 
zeta with polarizations functions (DZP) in combination 
with Troullier-Martins pseudopotentials. 51 Exchange and 
correlation are described by the PBE functional^ 2 which 
is close to PW91 used in the WFs calculations. The 
transmission functions obtained using the two different 
DFT-codes are seen to be in good agreement, especially 
in the important region near the Fermi-level. We have ob- 
served similar good agreements for the other Pt-CO-Pt 
structures considered in this work. This indicates that 
our results are independent of the basis set and other 
technical details related to the calculation of the self- 
consistent KS-Hamiltonian. 

In both cases (A and B) there is a pronounced reso- 
nance at the Fermi level, which clearly is responsible for 
the conductance of 1.5Go and 0.5Go, respectively. In 
fact this resonance is present for all the configurations 
investigated in this study, and it therefore represents a 
very general and robust feature of the transport through 
a Pt-CO-Pt contact. As the following analysis will show, 
the resonance is not due to the CO molecule alone, but 
rather is a (local) band structure effect related to the 
d-orbitals at the Pt apex atoms. 

To address the origin of the resonance, we perform an 
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FIG. 6: Typical transmission function for the upright bridge 
(configuration A in Fig. ^} and the tilted bridge (configura- 
tion B). The dotted line is the transmission function obtained 
using Siesta-DZP/PBE. Both transmission functions (A and 
B) have a resonance at the Fermi level which can be related 
to the density of states of the so-called group orbitals. 



analysis of the local electronic structure in the contact 
region by diagonalizing the Hamiltonian within the sub- 
space spanned by the WFs located at the CO molecule. 
The orbitals and eigenvalues obtained in this way repre- 
sent renormalized energy levels of the CO molecule in- 
cluding the effect of the coupling to the Pt leads. For 
all the considered contact geometries we find seven WFs 
located at the CO molecule, leading to seven renormal- 
ized molecular energy levels. For simplicity we focus 
on the tilted CO bridge in the following (configuration 
B), however, the conclusions are also valid for the up- 
right bridge configuration. Since CO has 10 valence elec- 
trons, the seven renormalized CO orbitals represent the 
five occupied and the two lowest unoccupied molecu- 
lar orbitals. The latter are the 2ir* orbitals which are 
known to be important for the chemisorption properties 
of C04L&£i By repeating the conductance calculations 
with the renormalized 2ir* orbitals removed from the ba- 
sis set, we find that the resonance at the Fermi level 
is completely gone and the conductance is reduced to 
~ 0.05Go- This allows us to focus exclusively on the 2ir* 
CO orbitals when analyzing the transport properties of 
the Pt-CO-Pt contacts. In the following we will refer to 
the 2ir* states as \a) and \b). The on-site energies of these 
renormalized orbitals are e a = 1.5 eV and Eb — 1.6 eV, 
respectively. The splitting of the levels is induced by the 
different couplings to the electrodes. 

Each of the molecular orbitals (MO) \a) and \b) 
give rise to one transmission channel through the CO 
molecule. If we neglect tunneling due to direct cou- 
pling between the Pt apex atoms, and neglect interfer- 
ence between the two transport channels, we can ana- 
lyze the problem by considering the transport through 
each MO separately. We do this using the well known 
single-level model for resonant transport. A particularly 
simple form of the single level model, is obtained by intro- 
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ducing the so-called group orbitals. 5 - The group-orbital 
for MO | a) (or \b)) of lead a (a = L,R) is defined as 
\da) — PaH\a), where H is the Hamiltonian, P a is 
the orthogonal projection onto lead a, and 1/V a , a a nor- 
malization constant. The name "group orbital" refers to 
the fact that \g%) consists of the group of states in the 
lead to which \a) is most strongly coupled. In fact, the 
coupling between \a) and any state in lead a orthogonal 
to \g£) is zero, while the coupling to the group orbital is 
given by V a , a — ( a \H\9a)- In the limit of strong asym- 
metric coupling, when << V£ , the transmission 
function takes the form 

T(s)=47r 2 Vl a p% a (e) Pa (e), (2) 

where p a (s) is the projected density of states (PDOS) of 
the MO \a) and p\ a (e) is the PDOS of the group orbital 
of the right lead in the absence of coupling to \a), i.e. 
calculated with Vr )( j = 0. The limit of strong asymmetric 
coupling is relevant for the tilted bridge configuration 
where V£/V£ « 0.1 for both MOs \a) and \b). The large 
asymmetry in the coupling strengths indicates that the 
Pt-C bond is much stronger than the Pt-0 bond. To 
illustrate the situation, the inset of Fig. [3 shows an iso- 
surface plot of |o) (transparent) together with its left and 
right group orbitals (solid). The latter consists mainly 
of d-like orbitals centered on the apex Pt atoms. The 
coupling strengths Vh,a an d Vr i0 are indicated. 




E-E F (eV) 

FIG. 7: (Color online) The PDOS of the 2tt* states \a) and |6> 
together with the PDOS of the group orbital in the weakly 
coupled right lead. The inset shows an iso-surface plot of 
the MO | a) (transparent) and its corresponding left and right 
group orbitals (solid). The PDOS for the MOs (circles and 
stars) is quite flat around the Fermi level, while the the PDOS 
of two group orbitals (full and dashed lines) both have a peak 
at the Fermi level. It is this peak that gives rise to the reso- 
nance in the transmission function. 

According to Eq. the energy dependence of the 
transmission function is determined by the product of 
the PDOS of the MO and the PDOS of the group or- 
bital in the weakly coupled lead. In Fig. we show the 
calculated PDOS for the MOs \a) and \b) together with 



the PDOS of the corresponding group orbitals \g^) and 
\g b ji). It is now clear that the transmission resonance at 
the Fermi level results from a corresponding peak in the 
PDOS of the group orbitals of the right lead, or, equiva- 
lently, from a peak in the PDOS of the d-states at the Pt 
apex atoms. The bare energies of the MOs at e a = 1.5eV 
and Eb — 1.6eV, respectively, are shifted upwards by the 
coupling to the Pt d-band, and can be seen as broad 
peaks in the PDOS at ~ 2.2eV. These peaks are also 
clearly visible in the transmission function in Fig. El 

The transmission resonance at Ep is thus caused by 
the properties of the isolated Pt lead, while the role of 
the 2ir* CO orbitals is to provide a flat background at the 
Fermi level and a peaked structure at ~ 2.2eV. Since the 
group orbitals for different configurations are very similar 
(always d-like) the determining factors in the transmis- 
sion functions stay almost constant, and this explains 
the robustness of the main features in the transmission 
function. 

It is well known, that the calculated HOMO-LUMO 
gap of CO is somewhat sensitive to the applied exchange- 
correlation functional4i However, since the transport 
properties of the investigated Pt-CO-Pt bridge involves 
only the tails of the PDOS of the CO orbitals and is 
dominated by the PDOS of the Pt leads, the exact po- 
sitions of the CO energy levels are not expected to be 
crucial. An accurate description of Pt is, however, im- 
portant and as discussed in Sec. lIIII our results for Pt are 
in good agreement with experiments. 



V. SUMMARY 

We have performed calculations within the framework 
of density functional theory for the mechanical and elec- 
trical properties of pure Pt nano-contacts and Pt contacts 
with a single CO molecule. 

For the pure Pt contacts, we obtain conductances 
traces which are in good agreement with experiments 
as well as other recent theoretical calculations. Our re- 
sults show that Pt point contacts have a conductance in 
the range (2.0 — 2.3)Go while that of short Pt chains is 
(1.3 — 2.0)Go- This provides a theoretical justification for 
assigning the peaks at ~ 2.1 Go and ~ 1.5Go in the con- 
ductance histogram for Pt to point contacts and chains, 
respectively. 

For the Pt-CO-Pt contact we identify an energetically 
stable configuration with the CO molecule providing a 
tilted bridge between two Pt apex atoms. Based on re- 
alistic DFT simulations of the creation of a conductance 
trace with the elastic response of the electrodes included 
through effective spring constants, we propose that the 
tilted CO bridge is responsible for the peak at 0.5Go ob- 
served recently in the conductance histogram for Pt-CO- 
Pt. We characterize and explain the main features of the 
transmission function for the Pt-CO-Pt contact in terms 
of the properties of the isolated CO molecule and the 
free Pt leads. The analysis shows that the conductance 
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to a large extent is determined by the local d-band at the 
Pt apex atoms and to a smaller extent by the 2ir* CO 
orbitals. 
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